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Abstract 

The peculiarities in the Bose-Einstein condensation of particles and quasiparticles 
are discussed. The difference between the condensation of conserved and unconserved 
^-j- ■ particles is analyzed. A classification of quasiparticles is given. The emphasis is made 

on the ability of particles and quasiparticles to condense. Illustrations include: general 
Bose-condensed atomic systems, such as ensembles of trapped atoms, Bose gases with 
conserved and unconserved number of atoms, vibrating atoms in double- well lattices, 
Holstein-Primakoff magnons, Schwinger bosons, slave bosons, and the condensation 
of singletons and triplons. The basic difference is that the system of particles, whose 
^ ■ total number is conserved, can form equilibrium as well as nonequilibrium condensates, 

while unconserved particles can condense only in a nonequilibrium system subject to 
external pumping supporting the density of these particles sufficient for their conden- 
sation. The examples of such a nonequilibrium condensation of unconserved particles 
are the Bose-Einstein condensation of excitons, polaritons, and photons. Elementary 
collective excitations, such as bogolons and phonons, being self-consistently defined, do 
not condense. Magnons cannot condense in equilibrium systems. Controversies, exist- 
ing in literature with regard to the Bose-Einstein condensation of some quasiparticles, 
are explained. Pushing a system out of equilibrium may favor the condensation of un- 
conserved quasiparticles, but suppresses the condensate fraction of conserved particles. 
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1 Introduction 



Bose-Einstein condensation has been a topic of intensive studies, both experimental and 
theoretical. There exists numerous literature on the problem, which can be found in the 
books [1-4] and review articles [5-18]. A variety of experiments have been performed with 
cold trapped atoms. In addition to atoms, whose total average number can be well controlled, 
hence, is conserved, one often considers Bose-Einstein condensation of unconserved particles, 
e.g., corresponding to excited bound states or collective excitations. The known examples are 
the Bose-Einstein condensation of excitons [19-23], polaritons [24-27], and, recently, photons 
[28]. 

Sometimes controversy arises, accompanying the discussion on the feasibility of realizing 
the Bose-Einstein condensation of unconserved particles. A very illustrative case of such a 
controversy is the discussion on the possibility of magnon condensation. In several works 
(e.g. [29-32]), the authors claimed the existence of magnon condensation in some equilibrium 
magnetic materials. Because of the known mapping [33,34] of spin systems to Bose gas, the 
existence of magnon condensation would mean the occurrence of Bose condensation for the 
collective excitations (bogolons) of Bose systems. Is it feasible that bogolons would condense, 
may be in a nonequilibrium Bose gas [35] in an external alternating potential? The problems 
of whether bogolons and magnons could condense are closely related. Mills [36] stressed that 
magnon condensation in equilibrium systems is impossible. Why this is impossible and could 
magnons condense in nonequilibrium systems? More generally, could elementary collective 
excitations become condensed and, if so, when? 

Understanding the necessary conditions, when particles and quasiparticles could form 
Bose-Einstein condensates, is a general problem that has become of great importance due 
to high current interest to Bose-condensed systems [1-18]. In the present paper, a careful 
analysis is given of such necessary conditions of condensation. The consideration is illustrated 
by the cases of the most common types of quasiparticles whose classification is also given. 
The suggested analysis resolves some controversies arising with regard to the possibility of 
Bose-Einstein condensation of quasiparticles. 

2 Particles versus quasiparticles 

Traditionally, one distinguishes elementary particles and quasiparticles. The distinction, 
however, can differ depending on the considered physical applications. 

In particle physics, an elementary particle, or fundamental particle, is a particle not 
known to have substructure; that is, it is not known to be made up of smaller particles. If 
an elementary particle truly has no substructure, then it is one of the basic building blocks 
of the universe from which all other particles are made. In the Standard Model, the quarks, 
leptons, and gauge bosons are elementary particles [37]. 

In statistical physics, one calls an elementary particle that whose structure is either 
absent or of no importance for the considered application. So, hadrons (mesons and baryons, 
such as the proton and neutron) and even whole atoms and molecules may be regarded as 
elementary particles, provided that their structure, if it exists, can be neglected for the 
studied phenomena. 

The notion of quasiparticles is not well defined. One often calls quasiparticles the dressed 
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particles whose properties are changed by the surrounding medium. For example, as an 
electron travels through a semiconductor, its motion is disturbed in a complex way by its 
interactions with all of the other electrons and nuclei. However it approximately behaves like 
an electron with a different mass traveling unperturbed through free space. This "electron" 
with a different mass is called an electron quasiparticle or Landau quasiparticle. Various 
single-particle and collective excitations are also termed quasiparticles. To make the consid- 
eration better defined, we give below a classification of different types of quasiparticles. 

2.1 Dressed particles 

Dressed particles are those that can be treated as elementary, but with the properties that are 
changed by surrounding, making them different from their properties in vacuum. Examples 
are: electron in a solid, electron or ion in plasma, proton or neutron inside a nucleus, polaron 
(electron with polarization of surrounding ions of a crystalline lattice), polariton (photon in a 
medium). Sometimes composite particles can mathematically be treated as elementary, such 
as atoms and molecules. A stable bipolaron (bound pair of two polarons) can be considered 
as one Bose particle. Strictly speaking, in quantum field theory, all particles interact with 
vacuum, being in that sense all dressed. The dressed particles can be bosons or fermions 
depending on their spin. 

2.2 Particle excitations 

Some quasiparticles behave as dressed particles, but differ from them by the necessity of 
supplying energy to the system for their creation. In that sense, such particle excitations are 
not stable and decay in an equilibrium system. For instance, a hole in the valence band of 
a semiconductor is a lack of an electron that has been excited to the conductance band. A 
hole is a fermion quasiparticle. Another example is an exciton in a semiconductor, that is a 
bound state of an electron and a hole, requiring for its creation photon excitation. Excitons 
can be treated as bosons. The number of excitations is not conserved, depending on the 
intensity of the applied external pumping. 

2.3 Collective excitations 

Collective excitations describe quantized fluctuations of an order parameter. Examples are: 
phonon (quant of mass density fluctuations), plasmon (quant of charge density fluctua- 
tions), magnon (quant of spin fluctuations), bogolon (quant of density fluctuations in a 
Bose-condensed system). Such collective excitations are usually described by linear equa- 
tions characterizing small deviations from a stationary state, because of which they are 
termed elementary collective excitations. All these collective excitations are bosons. They 
are not conserved. 

2.4 Nonlinear waves 

Another type of excitations can also be formed by the system as a whole, being in that 
sense also collective, but, contrary to the elementary collective excitations, being described 
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by nonlinear equations corresponding to large deviations from a stationary state. In clas- 
sical and quasiclassical picture such excitations are named nonlinear waves. As examples, 
we can mention solitons [38,39], bions (two bound solitons) [40,41], triple solitons (three 
bound solitons) [42], quantum vortices [17], and topological coherent modes [43-54]. Being 
quantized, such nonlinear excitations can be interpreted as quasiparticles. Usually they are 
not stable, except some one-dimensional solitons, arise only in nonequilibrium systems, and 
are not conserved. 



2.5 Auxiliary quasiparticles 

In quantum theory, particles and quasiparticles are represented by field operators of cre- 
ation and annihilation. The operators can be transformed resulting in the appearance of 
new operators that can correspond to other physical quasiparticles. Thus, the field operator 
of uncondensed particles in a Bose-condensed system can be decomposed into the linear 
combination of operators corresponding to collective excitations, bogolons. But it is not 
necessary that any operator transformation would result in new operators corresponding 
to some physical particles or quasiparticles. A transformation can be just a mathematical 
formula, convenient for calculations, but not necessarily related to new physical objects. 
Nevertheless, it is customary to associate each field operator, occurring in operator transfor- 
mations, with a quasiparticle, though this can be merely a fictitious particle being related to 
no real physical object. Such auxiliary fields and quasiparticles occur, e.g., in the Hubbard- 
Stratonovich transformation [55]. Other examples will be given below. Auxiliary particles 
can be either bosons or fermions, or even anyons, having intermediate fractional statistics 
[56]. 



3 Condensation of cold atoms 

There are two ways of treating Bose-Einstein condensation of particles whose total number is 
fixed. It is possible to consider a finite system of N particles in volume V. In a finite system, 
strictly defined phase transitions are absent. Following a kind of the Penrose- Onsager scheme 
[57], one needs to find the largest eigenvalue of the single-particle density matrix, which 
represents the number of condensed particles N Q . If the value n = N /N in thermodynamic 
limit remains finite, one says that Bose-Einstein condensation takes place. In that case, the 
order index of the single-particle density matrix becomes unity [58-60] . This way is extremely 
complicated and, for realistic models, can be realized only by means of powerful computers. 

Another way takes into account that, though strictly defined phase transitions, including 
Bose-Einstein condensation, occur only in thermodynamic limit, but, if the number of par- 
ticles in the considered system is sufficiently large, such that N 3> 1, then one can use the 
relations typical of infinite systems, where phase transitions are well defined. This essentially 
simplifies the consideration. In thermodynamic limit, Bose condensation is accompanied by 
spontaneous breaking of global gauge symmetry. Moreover, the symmetry breaking is nec- 
essary and sufficient for Bose-Einstein condensation [11]. In a finite system, the symmetry 
breaking is asymptotic, in the sense that it is approximate for finite N, and becomes exact 
for N — > oo. For a large number of particles N 1, the symmetry can be considered as 
broken. Then it is reasonable to take into account the symmetry breaking explicitly, which 
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makes the theory more transparent and calculations treatable. 

In the case of Bose condensation, the most convenient method of symmetry breaking is 
by means of the Bogolubov shift [61,62], when the Bose field operator is represented as the 
sum 

^(r,*)=» 7 (r,*)+^ 1 (r,*), (1) 

in which the first term is the condensate wave function and the second term is the field 
operator of uncondensed atoms. The latter two quantities are orthogonal to each other, 
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»7*(r,ftyi(r,f)dr = 0. (2) 



Uncondensed atoms are considered as normal, such that the statistical average of the operator 
of uncondensed atoms is zero: 

(^(r,t)} = 0. (3) 
Then the average of the total field operator (1), 

(j,(r,t))=r,(r,t), (4) 

gives the condensate wave function that can be treated as the order parameter. Its modulus 
squared yields the condensate density 

p (r,t) = \r ] (r,t)\ 2 , (5) 

which is normalized to the number of condensed atoms 

J h(r,t) \ 2 dr = N (t) . (6) 

Another normalization condition gives the number of uncondensed atoms 

(N^t)) = Nl (t) (7) 
that is the average of the number operator 

= J 4(r,t)M*,t) dr. (8) 

The sum of the numbers of condensed, (6), and uncondensed, (7), atoms gives the total 
number of atoms 

N (t) + N^t) = N . (9) 

Generally, the total number of atoms could depend on time, when there is the input of atoms 
or their loss. But in any case, this number is fixed at each moment of time by the given 
external conditions and cannot be created inside the system. 

Equation (3) can be rewritten in the standard form of a statistical condition by intro- 
ducing the operator 

A(t) = J [A(r,^I(r,t) + A*(r,t)Vi(r,t)] dr. (10) 



5 



Then Eq. (3) is equivalent to the condition 

<A(f)> = . (11) 

The above formulas are defined for an arbitrary system, equilibrium or nonequilibrium, 
provided that the condensate does exist, that is, when the condition of condensate existence 

Km Mi > (12) 

is valid. The system as such can be finite, but the thermodynamic limit in Eq. (12) serves 
as the check for the accurate definition of the condensate existence. This condition imposes 
constraints on the behavior of the distribution of uncondensed atoms and the spectrum of 
collective excitations [16,18] for an equilibrium system. The distribution of atoms n k over 
the quantum multi-indices k is to be such that, if the condensation is into the state labelled 
by k , then 

nk — > oo (k — > k ) . (13) 
This is equivalent to the requirement that the spectrum e k be gapless: 

e k ->• (k -» fc„) • (14) 

Another important point is that the cases, when the condensate exists or does not exist 
require different mathematics [63]. When the condensate exists, hence condition (12) holds 
true, then the action of all operators is defined on the Fock space ^(tpi) generated by the 
field operator ipi. But when there is no condensate, so that the condensate function r](r,t) 
and the condensate density (5) are zero, then the theory is defined on the Fock space J 7 ^) 
generated by the field operator ip, without any shift. These two spaces are orthogonal to 
each other [63]. It is incorrect to work in the space J 7 ^) and then make the Bogolubov shift 
(1) passing to the space J-{ip\). This is the standard mistake. As soon as there happens 
Bose-Einstein condensation, one has to work in the space J-tyi). 

Let H[r),if)i] be the energy Hamiltonian for the system. Because of the existence of three 
statistical conditions (6), (7), and (11), the grand Hamiltonian must contain three additional 
terms: 

H = H[ V , Vi] - IM>N (t) - fuMt) ~ A(0 • (15) 

Here /io and ji\ are the Lagrange multipliers guaranteeing the normalizations (6) and (7). 
And A(r, t) is the Lagrange multiplier cancelling the linear in ipi(r, t) terms in order to satisfy 
condition (11). The grand Hamiltonian (15) allows one to correctly define a representative 
statistical ensemble uniquely describing the Bose-condensed system [63-67]. 

The equations of motion can be represented by the variational derivatives of the grand 
Hamiltonian over r\ and ipi, which is equivalent to the Heisenberg equations of motion [18,68]. 
Thus, the equation for the condensate wave function reads as 

d . , / 5H \ 

In the case of an equilibrium system, this reduces to the equation 

5H 



. = . (17) 
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And, if the equilibrium system is uniform, the previous equation becomes 




= . 



(18) 



Evolution equations, together with the related initial conditions, define the functions 
N (t, /i , /ii) and N ± (t, /x , A*i)- Because of normalization (9), the multipliers /i and H\ are 
connected with each other and, generally, can depend on time. The condition of the gapless 
spectrum, for a nonequilibrium system, reads as 



These equations define the multipliers Ho{t,N) and fi\(t,N). 

The time-dependent solutions should tend to the corresponding stable stationary solu- 
tions N = N (fi , Hi) and N 1 = Aq(// , Hi) that, according to sum (9), satisfy the normaliza- 
tion N +Ni = N. The condition of condensate existence, in the form of the gap less-spectrum 
requirement (14), gives fii — Hi{no->N) by the Hugenholtz-Pines relation [69]. These four 
equations uniquely define all four quantities: N , JVi, fio, A*i through the parameters of the 
system in the stationary state. 

If one would assume the equivalence of the chemical potentials /i and this would 
overdefine the equations and would lead to inconsistences, as has been shown by Hohenberg 
and Martin [70]. Such inconsistences would lead either to the appearance of a gap in the 
spectrum or to the invalidity of the equation (16) for the condensate wave function. In turn, 
thermodynamic relations would break down and the poles of the single-particle and two- 
particle Green functions would be different, contradicting the results of Bogolubov [62] and 
Gavoret and Noziere [71], according to which the poles of these Green functions coincide, 
when the global gauge symmetry is broken. 

The Bogolubov canonical transformation 



introduces the Bose operators bk(t) of collective excitations, bogolons, where k is a quantum 
multi-index. This transformation makes it possible to diagonalize the grand Hamiltonian in 
the Hartree-Fock-Bogolubov approximation. The coefficient functions u k and v k are found 
from this diagonalization which yields the Bogolubov equations that also give the Bogolubov 
spectrum The bogolon number is not conserved and bogolons cannot condense. This 
is evident from relation (20). Really, if the bogolons would condense, resulting in nonzero 
(b k (t)), then the average {ipi(r,t)) would also become nonzero, which contradicts the basic 
definition (3). 

Thus, atoms, whose total number N is given, can form an equilibrium, as well as nonequi- 
librium, Bose-Einstein condensate. But the collective excitations, bogolons, whose number 
is not conserved, cannot be condensed. If in a nonequilibrium system one would notice 
the tendency of bogolons to condense, this would simply imply the necessity of redefining 
the condensate function of atoms in such a way that to preserve the principal definition 
(4), hence, avoiding the bogolon condensation. The same concerns other collective exci- 
tations that could arise in Bose-condensed systems, even if these excitations would have 
single-particle form [72]. The main point is that the collective excitations correspond to 



e k {t) -> 



(k -» fc ) . 



(19) 




(20) 



k 



7 



unconserved entities, so that they cannot condense. As is clear from the present section, bo- 
golons, that are elementary collective excitations in a Bose-condensed system, by definition 
cannot condense whether in equilibrium or in any nonequilibrium case. 



4 Conserved versus unconserved particles 

As an explicit illustration demonstrating the basic difference between conserved and un- 
conserved quasiparticles and the impossibility of equilibrium condensation for unconserved 
quasiparticles, let us consider boson quasiparticles with the effective spectrum 

oo k = Ak n , (21) 

where A and n are positive parameters and k is the modulus of a d - dimensional wave vector. 
This type of spectrum occurs in some problems of solid-state physics and some cosmological 
problems, where the matter with such a spectrum is named absolute stiff matter [73-76]. 
Quasiparticles are assumed to be in a large box of volume V, so that the sum over d- 
dimensional momenta can be replaced by the integral 




2k d ~ l dk 
(47r) d / 2 r(d/2) ' 



The uniform density is denoted by p = N/V . 



4.1 Equilibrium condensation of conserved quasiparticles 

Here we consider conserved quasiparticles, whose total number N is fixed, therefore the 
density is constant, p = const. The critical temperature, where the condensation starts, is 
found in the usual way [18], yielding 



T c = A 



(An) d l 2 T(d/2)np 
2T(d/n)((d/n) _ 



n/d 



(22) 



Taking into account that the gamma-function T(x) > for x > and the Riemann zeta- 
function ((x) < in the interval < x < 1 tells us that there is no condensation for d < n. 
When d — n, then T c = 0. And T c is positive only for d > n. In the latter case, for the 
condensate fraction below T c , one has 

jy / T \ d l n 



However, the formal existence of the critical temperature does not yet mean that the con- 
densate can really occur. In order to exist, it must represent a stable statistical system. For 
this purpose, the system compressibility is to be finite, which is equivalent to the finiteness 
of the stability ratio 

var(iV) . . 

< < oo (24) 
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for all N, including thermodynamic limit [7,16,18,77]. Here 

var(iV) = (N 2 ) - (N) 2 , N = N + Ni . 
Calculating ratio (24), we use the fact that 

var(iV) = var(iV 1 ). 

Then, for n < d/2, we get 

var(iV) _ t{d/n-l) ( T_\ d/n 



N 



C(d/n) \Tj 



(d > 2n) . 



(25) 



In the interval d/2 < n < d, similarly to Ref. [78], we employ the generalized Bose functions, 
in whose definition the integrals are limited from below by the minimal wave number k m i n = 
{p/NY' d . This gives 



var(iV) 2{d - n)N { ? n - d V d 



N 



(47r) d / 2 n(2n - d)T(d/2) \Ap n / d 



T 



[n < d < 2n) . 



(26) 



In particular cases, remembering that T c > exists only for d > n, we have the following. 
When d = 3 and n — 2, the critical temperature (22) is 



T r = An A 



P 



am 



2/3 



(d = 3, ra = 2) . 



(27) 



But ratio (26) diverges for large N, 



var( 



(28) 



hence, the stability condition (24) is not valid. 

When d — 2 and n — 1, the critical temperature is 



T c = 2A\ — 

V 7T 



(d = 2, n= 1) 



(29) 



But the stability ratio is strongly divergent, hence the system is unstable. 
When d = 3 and n — 1, the critical temperature (22) becomes 



T P = A 



TT 2 p 

03) 



1/3 



(d — 3, n — 1) , 



(30) 



where we have used the value ((2) = tt 2 /6. The stability ratio (25) is finite: 

var(iV) T 3 



N 



6pA 5 



(d — 3, n — 1) , 



(31) 



telling us that the system is stable. 

In this way, the uniform d- dimensional Bose gas of conserved quasiparticles with spectrum 
(21) can be Bose-condensed, provided that d > 2n. 
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4.2 No equilibrium condensation of unconserved quasiparticles 



If the total number of quasiparticles is not fixed, but they can be created, then, when 
decreasing temperature, no condensation occurs, but instead, the number of quasiparticles 
diminishes with temperature, which implies their decreasing density 

p _ 2T(d/n)C(d/n) /T\ d < n (32) 



(4n) d / 2 nT(d/2) \A 

This density is positive for d > n. In the particular cases, corresponding to those studied 
above, we have 

P = cff)(^) 3/2 (d = 3,n = 2), (33) 



P=l(L) 2 (d = 2,n=l), (34) 



(d = 3, n = 1) . (35) 



3 \2AJ 
C(3) (T 



7T 2 \A 



In all the cases, the density of unconserved quasiparticles diminishes with temperature, 
tending to zero. But there is no condensation. 

The similar situation happens in the case of equilibrium photons. In equilibrium, at 
lowering temperature, photons disappear instead of occupying the cavity ground state. But 
under an external pumping supporting photon total number inside an optical microcavity, 
photons can condense [28]. Such a photon condensation, creating an ensemble of coherent 
photons is, actually, analogous to the functioning of lasers. 

When the lifetime of created quasiparticles is sufficiently long, a quasiequilibrium con- 
densation of unconserved particles in a quasi-stationary state can be realized, as in the 
case of excitons [19-23,79]. But in absolute equilibrium, unconserved quasiparticles cannot 
condense. 



5 Vibrating optical lattices 

A very common quasiparticle, representing collective excitations, is phonon, characterizing 
density fluctuations caused by vibration of atoms in a lattice. In addition to the standard 
solid-state physics, different lattices can be created by means of laser beams forming standing 
waves where cold atoms are caught [15,16,80,81]. Phonons are unconserved quasiparticles. 
Hence, we should expect that they cannot condense in equilibrium. However, in some cases 
one talks on condensation of phonons. Here the analysis is given explaining that the phonon 
condensation in equilibrium systems does not occur. The formal appearance of phonon con- 
densation may happen in nonequilibrium situations, such as phase transitions. But defining 
phonons in a self-consistent way makes it possible to avoid phonon condensation in any 
system, equilibrium or not. 
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5.1 Double- well optical lattices 

For generality, let us consider the case of atoms that can happen in two states of a lattice. 
A straightforward example of such two-state atoms are the atoms in double-well lattices 
[82-92] and different double- well potentials [93-95]. Another case could correspond to two- 
band lattices [96]. The usual single-band lattice can be considered as a particular case of a 
two-band lattice. 

The derivation of the lattice model starts with the standard Hamiltonian 

H = J ^( v )H l (r)i)(Y) dr + 

+ X - J ^(r)^{r')®{r - r')^(r / )^(r) drdr' , (36) 
where <3>(r) is an atomic interaction potential and the single- atom Hamiltonian is 

ifx(r) = - ^ + V lat (r) + U ext (r, t) . (37) 

In the latter, Vi at (r) is a lattice potential and U ext (r,t) is an external potential, such as the 
trapping potential. Generally, the external potential can be time- dependent. 

Atoms, located at the spatial points r^, are enumerated by the index j = 1,2, ... ,N. 
Keeping in mind an insulating lattice, we can describe atomic states by the localized orbitals 
given by the Schrodinger equation 
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^n(r -Tj) = E nj ^ n (r - Tj) . (38) 
The orbitals form an orthonormal basis with the property 

C(r-r,)« r - r ,)* = W«. 
The atomic field operator can be expanded over the orbital basis, 

V>(r) = J2 c «^n(r - Tj) . (39) 



/ 



nj 



We assume that at each location there exists just one atom, which imposes the unipolarity 
condition 

^ ] C nj C nj = 1 5 c mj c nj = . (40) 

n 

The matrix element 

ET = J C(r - rO#i(r)V>n(r - r,) dr (41) 
of the single-atom Hamiltonian (37) is the sum 

E rnn = R rnn + yran + jjrnn (42) 
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of the terms 

K™ = J r m (r - r 4 ) (- ^ ^n(r - r,) rfr , 
VT = / - 'O^tolM' - r,) rfr , 

And the eigenvalues of Eq. (38) are 

E nj = K™ + V™ . (43) 
Then the matrix elements (42) can be written as 

E" 171 = 5 8- E ■ + U" 1 "- (44) 

We assume that only two lowest bands are occupied, so that n = 1,2. Then it is possible 
to introduce the pseudospin operators 



Sj - - (c^cij - c\ j c 2 j S j , - - (c\ j c 2 j - 4j c l^ 

S i = \ { c h c V + c lcy) , (45) 



which gives the atomic operators 

c ij c ij = g + S i ' c 2i C2 i = 2 ~ S i 



,t „ - 1 _L J 1 



c t . C2j .S; iS* , cj . Clj = 5? + iS'l . (46) 

The pseudospin operators satisfy the standard spin algebra. And the atomic operators can 
satisfy either boson or fermion commutation relations. 

To better understand the physical meaning of the pseudospin operators (45), we can 
introduce the left and right location operators 

l / N l 

-j= {Cij + C 2j ) , C jR - —j- 



jL = -7= (Clj + C2j) , C jR = -= (Cy - C 2 j) (47) 



describing the location of atoms in the left or right wells of a double well that is centered at 
the lattice point rj. Then the pseudospin operators (45) become 

Sj = \ (c] L c jR + c) R Cj^j , S i=~\ ( c )lCjR ~ c] R c jL ) , 

S i=\ { C W - c) R c jR ) . (48) 

This shows that Sj corresponds to tunneling transitions between the left and right locations, 
Sj characterizes the Josephson current between these locations, and Sj describes the pop- 
ulation imbalance between the wells. Thus, the operator Sj can be called the deformation 
operator, since, when atoms move to one of the wells, the spatial system symmetry is broken 
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and the system becomes deformed. For a finite system, such a deformation leads to the 
change of its shape. 

Substituting expansion (39) into Hamiltonian (36), we denote the matrix elements of the 
interaction potential as Afcj), Bfcj), and C(rjj), whose detailed definition can be found in 
Refs. [89,90,97] and where 



Also, we introduce the following notation: 



i 

In this way, Hamiltonian (36) transforms into 

H = NE + J2{^ - ^J + W) + 



21 



(49) 



(50) 



Note that the expression hj, playing the role of an external field, is nonzero only when the 
external potential U ext (r,t) is not uniform. As soon as U ext is uniform in space, the effective 



field hj is zero. 



5.2 Operators of atomic deviations 

Taking into account vibrational degrees of freedom can be done by introducing the operators 
of atomic deviations Uj by the equation 



r j ~ a i + u i > 



(51) 



where a, is a fixed vector related to a j-th lattice site. Since the interaction terms in Eq. 
(50) depend on the difference r^, it is convenient to use the relative deviation 



u 



Ui - u, 



(52) 



Treating Uj as small deviations around a,, we expand the interaction terms in Hamilto- 
nian (50) in powers of u^. In this expansion, we use the notations 



OA- 

Aij = A(a,ij) , Afj = 



da? ' 



_ d 2 Ajj 

lj ~ dafdcfi ' 



where 



I a «i I 



&4j — &i &j 



Analogous notations are employed in the expansions of other interaction terms. 



13 



Since the expression 



A ~ jv E Ai i ~ E Ai 



does not depend on the atomic indices, we have 



OA 



a/3 



d 2 A 



o, 



E^ = 2 E< E^ = °- 

U » j 

In the above summations, we keep in mind that % 7^ j, though this is not shown explicitly 
for the brevity of notation. Alternatively, we may set An = 0, In = 0, B u = 0. 
Following this way [97,98], we come to the Hamiltonian 



+EE 

ij a/3 



\ <<< + I B 



E« 



a u a - 
13 13 



^2 B if u ij u ij j SiSj- 

a/3 / 



a0 



(53) 



Vibrational and atomic degrees of freedom in the four-operator terms can be decoupled 
by the relation 



,a n fi \ Q-f c7 



fi\ I Qfl C7\ 



(54) 



while keeping untouched the lower-order terms. We suppose that the pseudospin correlation 
functions (S"S") weakly depend on the lattice indices, so that 



d 2 



dafdaf 

1 1 3{¥=i) 



E J ?^. 



31 ~ 



d 2 



3(&) 



da? da- 



3(&) 



Thus we come to the Hamiltonian 

H = EqN + H p h + H at + H int ■ 

Here, in the first term 

£o = £o + ! + ^EE w> - «^>) « 

ij a/3 



(55) 



(56) 



14 



The second term is the effective Hamiltonian of dressed phonons 



2m 2 

3 V af: 



with the renormalized dynamic matrix 



< = Aff - 2lg(S?Sfi + 2Bg(SfS?) . (58) 



The third term is the Hamiltonian of dressed atoms 



//»/ = E - q j s j) + E - 'a 13 ' 13 !!) > ( 59 ) 

with the renormalized atomic interactions 

= % + E ^ ««f - «X> > hi = h + E l t - «i «f > • ( 60 ) 

a/3 a/3 

And the fourth term describes atom-phonon interactions, 

^ = - 2 EEw, (6i) 

i a 

in which the effective force acting on an atom is 

i? = X)iW, (62) 
j 

where 

/•^ = I*StS] - B*S?S? . (63) 
5.3 Pseudophonon and phonon operators 

In order to introduce phonon operators, let us define in the standard way the phonon fre- 
quency Uks by the eigenproblem 

^EE< eikaiJ 4 = ^, (64) 

3 P 

where k is momentum, s = 1, 2, 3 is polarization index, and e^ s is a polarization vector. 
If one defines the phonon operators dk s by the usual relations 



'2N , 

ks 



E V mw fa e ks (d ks - e' 



ka, 



«^^E 7 =(*. + ^.)^. («) 
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the total Hamiltonian (55) does not become diagonal with respect to these operators d ks . 
That is, such operators are not true phonon operators, but rather some pseudophonon op- 
erators. To diagonalize the Hamiltonian, we need to accomplish an additional canonical 
transformation 

d ks = b ks + 2 - Y V F«e a ks e~^ . (66) 

V2N u ks ^mu ks 1^ 

This is equivalent [97] to introducing the phonon operators b ks by the nonuniform relations 



% 



Pi = " 



f2N 



Y V muj k s e ks (b ks - b!_ k ^j e tk a J 



where the additional vector term has the components 

A ? = 2 ££7?if, (68) 

3 P 

with the matrix 

a/3 _ J_ ST^ °ks c ks ik-SHj ( m \ 
ks Ks 

As a result of this nonuniform transformation (67), Hamiltonian (55) takes the form 

H = NE + H ph + H at + H 4 , (70) 

with the diagonal phonon part 

Hph = Y Uks \ h ls h ks + ^) , (71) 

ks ^ ' 

and with an additional four-atom interaction term 

V fg a/3 

Since the Hamiltonian is diagonal in operators b ks , it is these operators that should be 
accepted as real phonon operators. 



5.4 Shape distortion operator 

Operator (68), entering transformation (67), characterizes the shift of the deviation operator 
Uj caused by the existence of atomic degrees of freedom described by the pseudospin oper- 
ators S". Such a shift of atoms from their initial positions a, would lead to the distortion 
of the considered sample, because of which the operator A^ can be called the distortion 
operator. In the case of a finite system, this deformation results in the change of the sample 
shape. To illustrate more explicitly the form of the distortion operator, one needs to specify 
the phonon frequency. 
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If we take the phonon frequency in the Einstein approximation 

u ks = u , (73) 

then matrix (69) is 



This gives the distortion operator 



A" = , F t a . (74) 

Another common case is the Debye approximation for the frequency 

u ks = Q(k D - k)ck , (75) 
in which Q(k) is a unit step function and 

^=(67T 2 p) 1/3 

is the Debye wave vector modulus. Then for matrix (69), we get 

a/3 _ Si(fcpajj) 

13 2n 2 pmc 2 dij 

where Si (re) is the sine integral 
The latter possesses the asymptotic properties 



x sin t , 
dt . 



x 3 



Si(x) ~ x - — (x 0) , 

7T COSX 

bi(x) ~ — — (x oo) . 

Zi Ob 



The distortion operator becomes 



A: 



5.5 No equilibrium phonon condensation 

Real phonons are not conserved quasiparticles and, as is seen from the phonon Hamiltonian 
(71), 

(b ks ) = , (77) 

that is, equilibrium phonons cannot condense. At the same time, as follows from Eq. (67), 
the average deviation can be nonzero, since 

<«i> = (A?) • (78) 
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To explicitly show how the deviation can arise, let us resort to the mean-field approxi- 
mation for the correlation function 

{S?S?) = {S?){S?) (i^j). (79) 
Denoting the mean pseudospin 

CJ = {Sf), (80) 

we have for the average force (62) 

= E i^ c ' c i - f^ftCj) ■ (si) 

Therefore, the observed distortion (68) is, generally, nonzero. For instance, taking the 
phonon frequency in the Einstein approximation yields 

> = ^ftc] - ir;/ -:(-) . (82) 

3 

When C° is zero itself or does not depend on the index j, then distortion (82) is zero. But for 
C" 7^ and depending on j, there appears a nonzero distortion. For double wells, the average 
C? plays the role of an order parameter distinguishing a disordered phase with CJ = from an 
ordered phase where CJ ^ 0. The ordered phase exists at low temperatures below a critical 
temperature [89,90]. The ordering means the appearance of atomic imbalance between the 
wells of a double well. The order- disorder phase transition occurs if there is no external 
force hj defined in Eq. (49). When the latter is absent, the system Hamiltonian (55) is 
invariant under the inversion — > —Sj, which implies the existence of a degenerate ground 
state. This situation is analogous to the Jahn- Teller effect [99,100], where the existence of a 
degenerate ground state of a molecule leads to the distortion of the latter. 

The appearance of a nonzero distortion (82) could be associated with the condensation of 
pseudophonons described by the operators d^ s . However, it is evident that this condensation 
simply means the shift of the average location of each atom occurring at the point of the 
order-disorder phase transition. At this point the system is unstable. After the phase tran- 
sition has happened, one has to change the definition of atomic locations and, respectively, 
the definition of the deviation from this location. 

Alternatively, it is possible to define the atomic locations as the averages 

a, = fo) . (83) 
Then, in view of relation (51), one always has zero deviation 

(u,) = , (84) 
and, because of equality (78), zero distortion 

(A?) = • (85) 

In this case, the occurrence of the order-disorder phase transition is exhibited in the sharp 
variation, e.g., with temperature, of the average atomic location (83). But there is neither 
phonon nor pseudophonon condensation. 
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Condition (84) can be preserved in any given approximation by defining the grand Hamil- 
tonian 

H = H-'£ i X j -u j , 

j 

in which the Lagrange multipliers Xj are chosen so that to cancel the linear in the deviations 
uj terms, thus, guaranteeing the statistical condition (84). 

Even more, definition (83) can be extended to nonequilibrium systems by introducing 
atomic deviations through the equation 

Tj(t) =a,(0 + u j (0 , 

in which the first term is the average 

a,-(t) = ( rj (t)) . 

I that case, by definition, 

(u,-(f)> = , 

hence there is no symmetry breaking: 

(h s (t)) = o . 

Defining in such a self-consistent way atomic deviations makes phonon condensation absent 
in any nonequilibrium system. 



6 No equilibrium magnon condensation 

Collective excitations, characterizing spin oscillations, or spin waves, are magnons. Their 
definition can be done by the Holstein-Primakoff transformation [101,102] for the ladder spin 
operators 

q± qx I • qy 

This transformation for the spin operator S,-, located at r,-, reads as 



SI = ^J2S - b)b 3 b 3 , SJ = b) ^/2S - b)b 3 , 

S? = S- b]bj , (86) 

where S is the spin value and b 3 are Bose operators representing magnons. The square root 
of operators is to be understood as Taylor series in powers of the magnon density operator 

hj = b)bj , 

which implies the assumption on the small deviation S — S*. The eigenvalues of the magnon- 
density operator hj are 0,1,2,..., 2S. In the lowest approximation, 

S+ ~V2Sb 3 , SJ ~ V2S b+ . (87) 
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The Holstein-Primakoff magnons, by definition, characterize small spin fluctuations around 
the ground state, for which the average spin is directed along the z— axis, 

(Sj) = Se z , (88) 

so that there are no transverse components: 

(Sf) = (S?) = . (89) 

The latter means that 

(b s ) = . (90) 

By construction, the Holstein-Primakoff representation (86) presupposes that magnons 
describe small fluctuations around the saturated magnetization directed along the z-axis 
[101]. The appearance of magnon condensation would mean the rotation of the magnetization 
because of the arising nonzero transverse magnetization. But then the total magnetization is 
not directed anymore along the z-axis, as is assumed in the Holstein-Primakoff representation 
(86). That is, the latter is not applicable, so that magnons are not well defined. 

If in calculations one formally meets magnon condensation in an equilibrium state, when 
(bj) becomes nonzero, this implies that at that point a phase transition occurs, with the 
rotation of the average magnetization. Then one has to redefine the system ground state, 
for which the Holstein-Primakoff transformation (86) would be valid, since representation 
(86) is not applicable if magnons would be condensed. The rotation of system magnetization 
does not mean magnon condensation, but means the change of the system ground state. 

Magnons are not conserved and in equilibrium they cannot condense, (bj) = 0, similarly 
to the impossibility of equilibrium phonon condensation or to the equilibrium condensation 
of other unconserved quasiparticles. Impossibility of the equilibrium magnon condensation 
has been stressed by Mills [36]. 

There are interpretations of some magnetic effects in nonequilibrium systems as being due 
to magnon condensation. These can be the systems that are subject to a sufficiently strong 
external alternating field, whose frequency plays the role of a magnon chemical potential 
determining the nonequilibrium density of magnons [103-108]. Or these can be the spin 
systems prepared in a strongly nonequilibrium initial state and relaxing to equilibrium by 
essentially nonlinear dynamics [109-120]. The magnon condensation is manifested as the 
emerging state of precessing spins, where all spins precess with the same frequency and 
phase. This nonequilibrium magnon condensation is accompanied by the arising off-diagonal 
order, when the averages 

(Sf) oc e ioJot 

are nonzero, where Uq is the Zeeman frequency. The nonequilibrium magnon condensation is 
an example of dynamic coherence, contrary to static coherence related to the Bose-Einstein 
condensation of conserved atoms. 

As is mentioned in the Introduction, there exists a number of articles claiming the oc- 
currence of magnon condensation in equilibrium stable magnets. Such claims are based on 
confusion, when either the Holstein-Primakoff representation is used outside of its region of 
validity or transformations are employed to some auxiliary quasiparticles that have nothing 
to do with magnons. One could talk on magnon condensation at magnetic phase transitions 
[121-123], which does not contradict the above conclusion on the impossibility of equilibrium 
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magnon condensation, since the point of a phase transition is that where the system loses 
its stability, hence, an equilibrium system does not exist. But the occurrence of magnon 
condensation in the whole range of thermodynamic parameters, in an equilibrium system, is 
principally impossible. 



7 Operators of auxiliary quasiparticles 

I many cases, one accomplishes operator canonical transformations aiming at simplifying the 
solution of a problem. In such transformations, the newly introduced operators correspond 
to auxiliary quiasiparticles that do not necessarily represent physical particles, but can be 
just mathematical tools, not related to any physical objects. In many cases the statistics 
of such auxiliary quasiparticles are even not uniquely defined, and the latter can be equally 
treated as bosons or as fermions. Depending on whether the auxiliary quasiparticles are 
conserved and prescribed to be bosons, they can formally exhibit condensation in equilibium 
systems. In this and in the following section, some examples of auxiliary quasiparticles are 
given. 



7.1 Jordan- Wigner transformation 

As an example of a trasformation, introducing auxiliary quasiparticle, let us mention the 
mapping of spin operators of spin 1/2, for a one-dimensional chain, onto the fermionic 
operators. After this mapping, one expresses the fermionic operators in terms of bosonic 
ones using Abelian bosonization [124,125]. 
One considers the spin operators 

Si = \*j , (91) 



where Uj is the vector Pauli matrix. The Jordan- Wigner [126] transformation 

■> 1 \ / i~ x , \ 1 

n=l J \ n=l / 

maps the ladder spin operators onto fermionic operators dj. The latter do not correspond to 
some really existing particles or quasiparticles, but they just allow one to solve the studied 
one-dimensional problem. 
One may notice that 

S+ST = a] aj . (93) 
Thence, the last of transformations (92) can be written in the form 

This is in agreement with the general properties of ^-spin operators 



Sj ' Sj — - (S^Sj + Sj'S i ) + S-Sj — SijSj , 
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Sf = S+St + (S*) 2 - sj = S(S + 1) . 
The inverse transformation reads as 

< = exp (^-iir S nS?j S+ , aj = exp £ S+S^j S- . (94) 

The total number of the auxiliary quasiparticles is not conserved, since 

In this case, the auxiliary fermionic operators do not represent physical particles. If in the 
process of some canonical transformation one would employ a mapping to auxiliary bosonic 
operators, one could talk on the possible condensation of the related auxiliary quasiparticles. 
Recall, however, that equilibrium condensation can happen only for conserved quasiparticles. 

7.2 Condensation of Schwinger bosons 

Auxiliary quasiparticles can be introduced, for example, by means of the transformation 

S j = \Yl a \oP<*V a iV ( 96 ) 

a/3 

mapping the spin operator Sj at the spatial location rj to quasiparticle operators a JQ , in 
which a = 1, 2 enumerates two admissible states, for instance, spin projections. The operator 

oW» = a «p e i (7 = x, y, z) (97) 

7 

is a linear combination of the elements a^ af3 of the Pauli matrices, where e 7 is a unit vector 
in the direction labelled by 7 = x, y, z. Thus, the components of vector (96) are 



5 7 = JXXVv^ ' 



2 

a/3 

Explicitly, this transformation takes the form 

S+ = a) x dj 2 , S~ = a) 2 a n , S j = \ ( a ]i a ji - ^2^2) ■ (98) 
The quasiparticle operators satisfy the constraint 

J2a] a a ja = 2S. (99) 

a 

It is important to emphasize that be either bosonic or fermionic operators, thus, 

representing any type of quasiparticles, whether bosons or fermions. Transformation (98) 
was, first, introduced by Bogolubov [62,127], with the operators a,j a representing electrons, 
that is, fermions. Schwinger [128] treated dj a as representing bosons, because of which the 
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latter are called the Schwinger bosons [129]. If the index a is interpreted as spin projection, 
then the Schwinger bosons are fictitious quasiparticles that cannot exist in reality, which 
follows from their unphysical nature, as far as they are bosons possessing spin one-half. But 
they can be real particles when the index a corresponds to the enumeration of some other 
components. Also, pairons [130], that are tightly bound pairs of particles, could be connected 
with Schwinger bosons [131]. 

In view of constraint (99), the total number of quasiparticles in conserved: 

EE^W) = 25iV - ( 10 °) 

j a 

Being treated as bosons, conserved quasiparticles can condense, which has been used in 
the interpretation of magnetic transitions [132,133], though interpreting the magnetic order- 
disorder phase transition as the Bose-Einstein condensation of Schwinger bosons can lead to 
incorrect first-order transition, instead of the correct second-order one [134]. 

Schwinger bosons are conserved quasiparticles and are principally different from the un- 
consented Holstein-Primakoff magnons [134]. This is why the former can condense in equi- 
librium systems, while the latter cannot. 



7.3 Auxiliary slave bosons 

In the theory of strongly correlated Fermi systems, one uses the slave-particle approach based 
on the factorization of the electron field operator into two operators, one of which is bosonic 
and another, fermionic. There are two alternatives that are termed slave boson and slave 
fermion representations. 

In the slave boson representation [135-137], the electron operator Cj CT , describing the 
annihilation of an electron at the point rj with spin a, is formally split into the product 

Cia = b\U , (101) 

where bi is a spinless bosonic operator and ji„ is a fermionic operator with spin a. The 
operators satisfy the no-double-occupancy constraint 

blk + J2 &f* = 1 • ( 102 ) 

The bosonic operator ft, is said to represent an auxiliary charged spinless quasiparticle, 
named holon, while the fermionic operator f ia , an auxiliary neutral spin one-half quasiparticle 

spinon. 

By this definition, it is clear that the holons and spinons are fictitious quasiparticles. 
According to constraint (102), the total number of these auxiliary quasiparticles is conserved: 

= N > ( 103 ) 

where N is the total number of electrons. The number of holons separately, strictly speaking, 
is not conserved. But it can be treated as quasiconserved, when the number of spinons 
is limited by a value smaller than N. If the number of spinons is fixed by the doping 
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concentration, as is usually done, then holons become conserved exactly. Therefore, it is 
admissible to talk about holon condensation [138] or holon pair condensation [139]. 

In the slave fermion representation [140-142], one formally splits the electron operator 
into the product 

Cur = ftbia , (104) 

with a fermionic operator ft corresponding to spinless charged fermionic holons and a 
bosonic operator bi a describing neutral spin one-half bosonic spinons. Again, the no-double- 
occupancy constraint is imposed: 

£>U* + £7i = l, (105) 

a 

so that the total number of holons and spinons is conserved, 

E (z< 6 U*> + (ftfdj = n , 

being normalized to the total number of electrons. Here, the spinons are rather unphysical 
quasiparticles, being bosons with spin one-half. 

Fictitious quasiparticles, introduced by the slave particle representations, are not real 
particles, being just auxiliary tools for calculations. In these calculations, one has to pre- 
serve the imposed constraints that limit the efficiency of the representations [143]. Under 
conditions, when the bosonic quasiparticles can be made conserved, they are allowed for 
formally experiencing condensation. 



8 Condensation of singletons and triplons 

Auxiliary quasiparticles are also introduced in the theory of quantum antiferromagnets, 
as suggested by Sachdev and Bhatt [144]. Generally, these quasiparticles can be defined 
either as bosons or as fermions. Treating them as bosons makes their condensation possible, 
provided they are conserved. 



8.1 Bond-operator representation 

Suppose there are two spins one-half, with the related operators Si and S 2 . This pair of 
spins can be treated as a dimer connected by bonds. Sachdev and Bhatt [144] suggested for 
the pair of the spins the bond-operator representation 



•>i — — J sh a + t^s i £ a /3 7 ^gt 7 J , 
1 V M J 



\ lsh a + tls + iJ2 

\ 01 



01 

t. 



S% = - 2 I » T *a + V a S + % ^ Sa^t^ ] , (106) 

expressing the spins through the quasiparticle operator s, characterizing the singlet state, 
and the operators t a , with a = x,y,z, describing the triplet state. Here, is the totally 
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antisymmetric unit tensor. The corresponding auxiliary quasiparticles can, therefore, be 
called singletons, and triplons, respectively. Explicitly, Eqs. (106) read as 

Si — — [s t x + t x s — % {t y t z — t z ty)\ , S\ = — [s ty + t y s — i (t z t x — , 



Si — — [sH z +t\s i [t x t y t\tx)\ , 



for the first spin, and as 



S2=-\ [Sh x +^,3 + 1 (t\t Z - t\ty)] , St = -\ [Shy +t\ S + l (t% ~ tit,)] , 



[Sh z +tls + i{tlty-tlt x )} , 



for the second spin. This representation could be considered as a generalized Schwinger 
representation (96) extended to the pair of spins. 

The restriction on the physical states to be either singlets or triplets leads to the constraint 



(107) 



with a — x,y, z. 

The spin operators satisfy the standard commutation relations 

Si, S*fj = j y ] £gP~/Si , S%, 5*2 = ^ y 1 £ a/3fS2 , S" , S^ = 



(108) 



and are normalized as 



Si = S 2 2 = S(S + 1) 



4 



(109) 



Under these conditions, the correct algebra of spin operators can be reproduced by either 
bosonic or fermionic commutation relations for the operators of the auxiliary quasiparticles. 
The statistics of these quasiparticles is not uniquely defined; they can be either bosons or 
fermions. If one interprets them as bosons, then the bosonic commutation relations are 
assumed: 



[s, S t] = i , 



t a , t 



S, 



a/3 , 



[S, t a ] = [s, tl] = 



The scalar product of two different spin operators is 

Si • S2 



V a 



S^S 



and the summary ^-component is 

Si + 5*2 = ~~ i {fyy ~~ tytx) ■ 
One introduces the following canonical transformation 

t X = (t + + t-) , ty = - (t+ - t-) 



tz — tn 



(110) 



(111) 



(112) 



(113) 
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The new operators satisfy the bosonic commutation relations 

fc, i] = 5, u , [t„ t v ]=0, (114) 

where /i, v — +, — , 

The constraint (107) takes the form 

s f s + ^4^ = l, (115) 

V 

with v = +, — ,0. The z-component sum (112) becomes 

Sl + Sl = tlt--tlt + . (116) 

8.2 Lattice of dimers 

A lattice, consisting of two spin sublattices corresponding to two components, can be treated 
as a lattice of dimers. Each pair of spins in a dimer can be described as is explained in the 
previous subsection. Constraint (115) takes place at each lattice site enumerated by the 
index j — 1, 2, . . . , TV, 

V 

where v — +, — , 0. This implies that the total number of quasiparticles is conserved: 

e(( s K)+E^W)) = n - ( 118 ) 

One usually assumes [144] that the singletons are completely condensed, so that Sj can 
be replaced by the nonoperator real quantity s equal to the average 

* = (*i), (H9) 

which does not depend on the site index because of the ideality of the lattice. 

Another standard assumption is that the main contribution among the triplons is due to 
only one type of triplons denoted by 

t^t J+ . (120) 
Under these assumptions, the bond-operator representation reduces to 

SI, = - \ t>fy = . (121) 
And the normalization condition (118) becomes 

Ns 2 + Y^ifyj) = N . (122) 

j 

Taking account of the lattice periodicity simplifies normalization (122) to 

s 2 + (t)t 3 ) = 1 . (123) 

The number of singletons is defined self-consistently by requiring the system stability. Then 
triplons are conserved quasiparticles that can experience Bose-Einstein condensation. 
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8.3 Order parameters 



There can exist two phases in the system, the normal phase of uncondensed triplons and that 
of condensed triplons [144]. The order parameter distinguishing these phases is the average 
(tj). The normal phase, where triplons are not condensed, corresponds to magnetically 
disordered phase, while the condensed triplon phase is magnetically ordered [144]. The 
second order parameter is the absolute value of magnetization 



M = 



N 



reduced to the number of lattice sites. With representation (112), we have 



M 



And representation (116) gives 



M 



1 

N 



1 

N 



-t ] t ■ ) 



t)_tj. 



Resorting to simplification (121) yields 

QX qy qy 

>=>lj - Z 2j , D]_j - Z 2j , 

Finally, using normalization (123), we get 



qz 



qz 



1-s 2 



-tit 



j" 3 ■ 



(124) 



(125) 



8.4 Disordered phase 

The system Hamiltonian was initially a functional of the spin operators and S 2 j. With the 
bond-operator representation (106), it became a functional of the quasiparticle operators Sj 
and tj a , with a — x,y, z. After the canonical transformation (113), the Hamiltonian became 
a functional of Sj and tj U , where v — +, — , 0. With the following simplifications, resulting in 
Eqs. (121), the Hamiltonian acquired the form of a functional of s and tj. The sequence of 
these transformations can be symbolically depicted as 

H[S lj: S 2j ] — > H[ Sj , t ja \ — )• H[ 8j , t jv \ — )• H[s, tj] . (126) 

In the phase, where triplons are not condensed, one has 

(tj) = . (127) 

To guarantee the normalization condition (122), one has to define the grand Hamiltonian 

H = H[s,tj] - /i (^Ns 2 + > (128) 
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with the Lagrange multiplier fj, being the system chemical potential. The number of single- 
tons is defined by minimizing the system thermodynamic potential, which leads to 




(129) 



This is a necessary condition for the system stability. 



8.5 Triplon condensation 

The Bose-Einstein condensation of triplons implies the appearance of the order parameter 

V = (tj) ^ , (130) 

which does not depend on the site index due to the lattice periodicity. In order to correctly 
describe the condensed phase, we introduce the Bogolubov shift 

tj = r, + Cj . (131) 

Thus, the Hamiltonian now is a functional of s,rj,Cj, which means the Hamiltonian trans- 
formation 

H[s,tj]^H[s, 77, cj] . (132) 
The following numbers of quasiparticles are defined: the number of singletons 

N s = s 2 = Ns<2 > ( 133 ) 
j 

the number of condensed triplons 

No = H 2 = Nr ? > ( 134 ) 

j 

and the number of uncondensed triplons 

= = (^i> » ( 135 ) 

where the number operator of uncondensed triplons is 

The total number of quasiparticles, according to Eq. (123), is conserved: 

N s + N + Ni = N . (136) 

To take into account these normalization conditions (133) to (135), it is necessary to 
introduce the grand Hamiltonian with the required number of Lagrange multipliers, which 
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defines a representative ensemble [16,18,145]. Following the general approach, as applied to 
Bose-condensed systems [64-68], we need to introduce the grand Hamiltonian 

H = H[s, rj, cj] - /2 S N S - fi N - . (137) 

The Lagrange multipliers /io, and /ii guarantee the validity of the normalization conditions 
(133), (134), and (135), playing the role of partial chemical potentials. The system chemical 
potential is defined by the variation of the system free energy F, which gives 

dF _ dF dN s dF dN OF dN x 

^ ~ M ~ 1m s 1m + m~ 7m + dN~ 1 lm ' (138 ^ 

This, with the partial derivatives 

dF dF dF . . 

» s = dN s > «> = w - o , ^i = ^. ( 139 ) 

yields the system chemical potential 

/i = fi s n s + fi n + ii x n x , (140) 

in which the quasiparticle fractions are defined as 

dN s dN dN, 

n *=div' no= 9iv' ni= 9iv- (141) 

Not all partial Lagrange multipliers /i s , /io, and jii are independent. Their relation follows 
from the condition of stability 5F = 0, which means 

dF rAr dF rAr dF rAr n 
5N S + — — 5N + — — 5N ± = . 



dN s dN u dN t 

In view of Eqs. (139), this is equivalent to 

fi s 5N s + fi 5N + firfNi = . (142) 
Because of normalization (136), one has 

5N S + 5N + 5N 1 ^0. (143) 

Using this, Eq. (142) reduces to 

(/i - ii s )SN + (/ii - Hs^Nx = 

and further to 

(/i - /i s )n + (/ii - (J, s )ni = . 

The latter gives 

/i s = /i ° n ° + /iini . (144) 
n + rii 

Substituting this fi s into Eq. (140) results in 

H = Us ■ (145) 
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The quasiparticle fractions 



n s = — = s, ™o = ^, n i = ^2^( c ) c i) ( 146 ) 

j 

satisfy the normalization 

n s + n + ni = l. (147) 
Hence, the fraction of condensed triplons is given by the equality 

n = 1 - m - s 2 . (148) 

Magnetization (125) is 

M = 1 - s 2 = n + n x . (149) 

The stability conditions, requiring the minimization of the thermodynamic potential, 
lead to the equations 

dH\ /dH\ . . 

0, ( — )=0, (150) 



dn s I \ dn i 

where the multipliers /i s ,/io,/^i are fixed. The first of these equations is equivalent to Eq. 
(129). Finally, the chemical potential fi± is defined by the condition of the triplon condensate 
existence [16,18], which is equivalent to the requirement of a gapless spectrum, following from 
the Bogolubov [62] and Hugenholtz-Pines [69] theorems. 

In this way, all quantities are uniquely defined. It is worth stressing the necessity of 
taking account of anomalous averages (tjtj). Neglecting these averages is mathematically 
unjustified, distorts thermodynamic relations, and spoils the phase transition order, mak- 
ing the Bose-Einstein condensation a first-order transition, as has been explained in Refs. 
[16,18,67,146]. While the proper account of the anomalous averages preserves the correct 
second order of the phase transition for Bose atoms [16,18,147,148], as well as for triplons 
[149,150]. 

Another principal point is the necessity of employing the representative ensemble with 
the grand Hamiltonian (137) in order to uniquely define a stable Bose-condensed system. 
For this, one has to introduce in the Hamiltonian the terms with the appropriate Lagrange 
multipliers. Forgetting some of these terms makes the system unstable. The majority of 
articles on triplon condensation do not use the correct grand Hamiltonian, hence, consider 
unstable systems. 

Let us emphasize that condition (129) fixes the number of singletons. As far as the 
total number of all quasiparticles N has been fixed from the beginning, being just the 
number of lattice sites, the number of triplons is also fixed. Thence, triplons are conserved 
quasiparticles and can condense. Being conserved quasiparticles, triplons are principally 
different from unconserved magnons. The latter cannot experience Bose condensation in an 
equilibrium system. 

It is also important to remember that triplons are auxiliary quasiparticles that, strictly 
speaking, are not precisely defined, since in the bond-operator representation (106), they 
can be considered either as bosons or as fermions. Choosing them as bosons is an arbi- 
trary decision. They could equivalently be chosen as fermions. So, the introduction of 
bosonic triplons is just a mathematical trick convenient for calculations, but the physical 
nature of these triplons as bosons is rather ambiguous. Contrary to this, the definition of 
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magnons through the Holstein-Primakoff representation (87) uniquely prescribes to magnons 
the Bose-Einstein statistics. Magnons are bosons by their birth. And magnons describe the 
straightforward physical reality characterizing spin oscillations. 

Unfortunately, many authors confuse triplons and magnons and write about magnon 
condensation in equilibrium dimer lattices. Of course, such a terminology, as has been 
stressed by Mills [36], is basically wrong. Triplons have nothing to do with magnons. The 
former are conserved quasiparticles and can condense in equilibrium, while the latter are not 
conserved and can never display equilibrium condensation. 



Unconserved quasiparticles cannot condense in equilibrium, although some of them can con- 
dense in nonequilibrium systems, e.g., excitons, polaritons, and photons. Coherent motion of 
transverse spins can also be interpreted as nonequilibrium magnon condensation. However, 
one should keep in mind that the Holstein-Primakoff transformation is not valid for strongly 
nonequilibrium spin systems, because of which it cannot not be used for the latter. 

Nonequilibrium is favorable for the possibility of condensation of unconserved quasipar- 
ticle. But, contrary to this, nonequilibrium perturbations usually destroy the condensate of 
conserved particles. In this regard, it is illustrative to study what happens with the Bose 
system, prepared in the Bose-condensed state, after the system is subject to the action of 
increasing external perturbation driving the system far from its equilibrium state. 

9.1 Coherent topological modes 

Trapped atoms provide a very convenient object for studying well controlled transformations 
of an equilibrium Bose-condensed system into a highly excited state. There are different 
ways of exciting a Bose-condensed system by external fields. All of these ways can be 
mathematically classified into two categories. One way is to act on the atomic density 
by applying an external field that is characterized by a time- dependent potential V(r,t). 
Another method is to make the atomic scattering length a s (t) a function of time by invoking 
the Feshbach resonance techniques. Generally, the scattering length can be made a function 
of time and space. These two methods lead to similar consequences. 

When the amplitude of the perturbing field is small, elementary collective excitations are 
produced, representing density fluctuations around the given equilibrium state. For short 
wavelengths A, much shorter that the effective system size L, collective excitations in a 
nonuniform system are described by the Bogolubov spectrum 



where c(r) is local sound velocity. For long waves, such that A ^> L, the spectrum of 
collective excitations is discrete and depends on the trap geometry [1-12]. These elementary 
collective excitations do not change the condensate fraction. 

But even rather weak perturbations may drive the system far from equilibrium, if the 
external field alternates with a frequency that is in resonance with a transition frequency 
between two energy levels of the trapped system. Recall that, for a trapped system, the 



9 Strongly nonequilibrium condensates 
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spectrum of the stationary Gross-Pitaevskii equation is discrete and, because of atomic 
interactions, it is not equidistant. This allows one to tune the frequency of the external 
alternating field close to one of the transition frequencies. 

Coherent topological modes are the solutions of the stationary condensate-function equa- 
tion [16,18,43,53,54]. At temperatures close to zero and asymptotically weak interactions, 
the condensate-function equation reduces to the stationary Gross-Pitaevskii equation 

#[y? n ]y? n (r) = E n ip n (r) , (151) 

with the Gross-Pitaevskii Hamiltonian 

H[<p] = - ^ + U(r) + iV $oM 2 , (152) 

in which U (r) is a trapping potential, N is the number of condensed atoms, and 

t , a s 

$ = 47T — 

m 

is interaction strength, with a s being scattering length. The modes are called coherent, since 
they correspond to the Bose-condensed part of the system, which, by definition, is coherent. 
And they are termed topological because the corresponding atomic density is topologically 
different from the density of the ground state. 

Suppose that, in addition to the trapping potential, the system is subject to the action 
of an alternating external potential 

V(r, t) = Vi(r) cos(wt) + V 2 (r) sm(ut) , (153) 

whose frequency is tuned close to a chosen transition frequency 

u 21 = E 2 -E 1 . (154) 

This implies that the resonance condition 



W 2 1 



< 1 (Aco =u- u 2 i) (155) 



is valid. The problem is analogous to that of resonance transitions of atoms, with the 
principal difference in its nonlinearity caused by atomic interactions. 

Dynamics of the system is described by the time-dependent condensate-function equation 
[16,18,53] or, at zero temperature and weak interactions, by the Gross-Pitaevskii equation 
that is nothing but a nonlinear Schrodinger equation 

i ^ v?(r, t) = {H[<p} + V(r, t)}<p(r, t) . (156) 

The solution to this equation can be represented in the form of the expansion over the 
coherent modes, 

V9(r,t) = ^ Cri (t)v9„(r)e-^. (157) 
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The fractional mode populations are given by the expressions 

p n (t) = |c n (t)| 2 . (158) 

Under the resonance condition (155), it is feasible to excite the coherent modes employing 
rather weak external fields. Similarly, by several alternating fields, one can generate several 
coherent modes [51,52]. If we increase the amplitude of the perturbing external fields, then 
the topological coherent modes, in addition to the direct resonance, with oj = oj 2 i, can also 
be generated by other processes, such as harmonic generation under condition 

nu = u>2i (n — 2, 3, . . .) , (159) 

parametric conversion, when 

ui ± u 2 = u 2 i , (160) 

and other combined resonances, when 

^2niUi = U2i (rii = ±1, ±2, . . .) . (161) 

i 

Thus, the increasing perturbation can generate a variety of different coherent modes. And, 
with a sufficiently strong perturbation, strict resonance conditions are not required. 

Among all these modes, there are quantized vortices. Their generation in a trap, per- 
turbed by external fields, can be done so as to produce vortices with a chosen circulation, or 
the pairs of vortices and antivortices can be generated. One often employs cylindrical traps 
characterized by a transverse frequency u±_ and longitudinal frequency u z , with the aspect 
ratio 

a^=/>V, (162) 



where the related oscillator lengths are 



I I = -=L= ■ 



The strength of atomic interactions can be described by the dimensionless coupling parameter 

g = 4nN . (163) 

The basic vortex, having unit circulation [4], possesses the energy 



R 



where the notation is used for the Thomas-Fermi radius squared 



1 fl5 ^ 2/5 



Rtf = j- a 9 



and the healing length 



v/2mp$ ' 
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The vortex energy can be written in the form 



9u± 

Uvor = ( ag y/ 5 ln (0.8«5f) . (165) 

The number of atoms iV in a trap is usually so large that the coupling parameter (163) 
is a large number g >> 1. Then the dependence of the coherent-mode transition frequencies 
on the coupling parameter is such that 



oc {agf* , 



hence, these energies increase with g. While the energy of the basic vortex (165) decreases 
with g. This is why the basic vortex is the most energetically stable among all coherent 
modes, so that rather intensive perturbation creates mainly such basic vortices [5,18]. The 
core of a vortex is the region with condensate depletion. Thence, vortex generation depletes 
the condensate fraction. 



9.2 Creation of separate vortices 

When the strength of external perturbations continues growing, this pumps into the system 
additional energy. The action of an alternating potential V(r, t) increases the system energy 
so that the injected energy per atom is 



E inj = -jy J p(r,t) 



dV(r,t) 



dt 



drdt . (166) 



Let the pumping potential have the amplitude A and frequency u. Then the injected energy 
(166), during the time interval [t,t'], approximately is 

E inj ~ Au(t - t') . (167) 

A single vortex is created when the injected energy is of the order of the vortex energy, 
Ei n j ~ u vor , which implies 

Au(t - t ) ~ OJvor ■ 

Therefore the crossover line of a single vortex creation, on the amplitude-time plane, is given 
by the expression 

^ W^) ■ (168) 

Several vortices, of the number N vor , can be generated when the injected energy reaches 
the value E in j ~ N vor u vor , hence when 

Au(t — ti) ~ N vor u vor . 

The more vortices are produced, the more condensate fraction is depleted, though the system 
remains yet superfluid. 
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9.3 Regime of quantum turbulence 

Producing more and more vortices (and antivortices), one can reach the state where the 
vortices form a random tangle. Such a random tangle of vortices signifies the appearance of 
the state of quantum turbulence [151-158]. The number of vortices, when the random tangle 
is being formed can be estimated [18] as 

N vor ~ | , (169) 

where the effective trap parameters are 

7 1 /j2 i \l/3 1 / o \ V3 

This defines the crossover line of the starting turbulence 

Atur ~ - l f vor - . (170) 

In the turbulent regime the condensate fraction is getting more and more depleted. But the 
system can yet be treated as superfluid. 

9.4 Spatial condensate granulation 

With the developed turbulence, the condensate fraction becomes small. And finally, the 
condensate breaks into spatially separated pieces, or grains, that are immersed into the 
normal phase without condensate. Then locally, in each grain, there is condensate. However, 
since the grains are spatially separated, superfluidity through the whole volume cannot arise. 
Generally, such a state is called heterophase, as far as it is formed by a mixture of different 
phases, randomly distributed in space [145,159-162]. In the present case, it can be termed 
granulated condensate [163]. The typical size of grains with condensate is given by the 
localization length 

= ^kji = te) k ■ (m) 

These grains are mesoscopic in size, in the sense that the above localization length is larger 
than the mean interatomic distance a, but shorter than the effective trap size: 

a < hoc < lo ■ 

The granulation starts when the injected energy is so large that the localization length 
becomes of the order of the trap size, 

hoc ~ k , Einj ~ w . (172) 
This defines the crossover line of granulation 
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Increasing the injected energy diminishes the grain size (171). And the condensate is 
completely destroyed, when the localization length (171) reduces to the order of the mean 
interatomic distance: 

hoc ~ a , E inj ~ wo J— . (174) 



The corresponding crossover line 



Anor ~ -J^-, (175) 

u{t -t 3 ) V a 



shows where the granulated condensate is destroyed and the atomic cloud becomes a normal 
system without condensate. Since this system is strongly nonequilibrium, it can be named 
chaotic fluid. 



9.5 Amplitude-time phase diagram 

The above discussion can be summarized with the phase diagram in Fig. 1 on the plane of 
the pumping field amplitude versus pumping time. The sequence of these states, except the 
chaotic fluid, has been observed in experiments [164-167]. The state of chaotic fluid has not 
yet been reached. The details for the phase diagram can be found in Refs. [166,167]. 

In experiments [164-167], trapped atoms of 8 7Rb have been treated, having mass m = 
1.445 x 1CT 22 g and scattering length a s = 0.577 x 1CT 6 cm. The cloud was cooled down to 
low temperatures, much lower that the critical temperature of condensation T c = 276 nK, 
so that almost all N = 2 x 10 5 atoms were condensed. A cylindrical trap was used with the 
trap characteristics 

u± = 2vr x 210 Hz = 1.32 x 10V\ co z = 2ix x 23 Hz = 1.45 x 10V\ cu = 0.63 x 10 3 s"\ 

l ± = 0.74 x 10~ 4 cm , l z = 2.25 x 10~ 4 cm , l = 1.08 x 10~ 4 cm . 

This is an elongated trap with the aspect ratio a = u z /uj_ = 0.11. 

The effective atomic volume V e ff, average density p, mean interatomic distance a, and 
gas parameter 7, respectively, are 

v *ff = nl ± 2l z = 27r/ o = °- 78 x 10~ n cm 3 , p ~ 2.55 x 10 15 cm" 3 , 

a = p- 1/3 = 0.73 x 10" 5 cm , 7 = p 1/3 a s = — = 0.079 . 

a 

The dimensionless coupling parameters are 

g = AtiN y- = 1.96 x 10 4 , ag = 2.16 x 10 3 , 

which implies strong effective atomic interactions. 

The trapped system of condensed atoms was subject to the action of a modulating field 
of frequency 

u = 2tt x 200 Hz = 1.26 x 10 3 s _1 , 

hence, with the modulation period T mo d = 2n/u = 5 x 10 -3 s. The action of this perturbation 
lasted during the modulation time t mod = 0.02 — 0.06 s. 
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The time of local equilibration is 

hoc = t = 0.93 x i(T 4 s . 

npa s 

Then the relation between the characteristic times is 

Woe T m od ^ t mo d . 

This means that the atomic system, during the process of modulation, has always been in 
the state of local equilibrium. 

The healing length, defining the vortex core, is 

f = ; 1 = 0.52 x l(T 5 cm . 
\/87rpa s 

Therefore the characteristic lengths are related as 

a s < £ ~ a < l . 

The vortex energy is u vor = 0.41 x 10 3 s -1 . The turbulent regime starts when the number 
of the generated vortices reaches the value 

N vor = — ~ 20 , 

which is in good agreement with experiment [166,167]. 

The consideration of the present section illustrates how the system of conserved atoms, 
initially prepared in an almost pure Bose-condensed state, and being subject to the pertur- 
bation by external alternating fields, goes through a sequence of nonequilibrium regimes, 
ending with a normal chaotic fluid without condensate. The initial condensate fraction has 
been almost 100%, while the final state contains no condensate. Nonequilibrium external 
perturbation depletes the condensate fraction of conserved particles. 

A challenging problem would be to investigate experimentally the opposite process of 
equilibration of an excited atomic system, from the state of chaotic fluid back to the Bose- 
condensed state. Such problems of equilibration of finite quantum systems attract now high 
attention, as can be inferred form the review articles [168-170]. 

10 Conclusion 

For the phenomenon of the Bose-Einstein condensation, the distinction between particles 
and quasiparticles is secondary. The most important is whether their number is conserved 
or not. Conserved particles and quasiparticles can experience Bose-Einstein condensation 
in equilibrium systems as well as can form nonequilibrium condensates. But unconserved 
quasiparticles cannot exhibit equilibrium condensation. However, they can display conden- 
sation in nonequilibrium states, when an external pumping creates a sufficient number of 
quasiparticles. Examples are nonequilibrium condensates of excitons, polaritons, and pho- 
tons. 
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It is possible to interpret the shape and Jahn- Teller phase transitions as condensation of 
phonons. However this interpretation is useful solely at the phase transition points, where the 
system is actually unstable. After the phase transition has occurred, there are no condensed 
phonons, but it is just necessary to redefine the mean atomic locations. It is possible to 
avoid talking on phonon condensation by interpreting shape transitions as a sharp variation 
of the mean atomic locations. Self-consistently defined phonons never condense, whether in 
equilibrium or nonequilibrium systems. 

Among quasiparticles, there are those that enjoy direct physical meaning, such as phonons 
and magnons. And there are auxiliary quasiparticles, appearing in the process of mathemat- 
ical transformations and having no precise physical meaning, for instance, Schwinger bosons, 
slave bosons, or singletons and triplons. These auxiliary quasiparticles even are not uniquely 
defined as bosons. If, by force, they are assumed to be bosons, they can display equilibrium 
condensation, provided they are conserved. 

When considering quasiparticle condensation, it is principally important to employ a 
correct terminology, not confusing conserved with unconserved quasiparticles. For example, 
triplons are auxiliary conserved quasiparticles that are formally allowed to condense in equi- 
librium. And triplons must not be confused with magnons which are physical unconserved 
quasiparticles that cannot experience equilibrium condensation. 

The appearance of average transverse magnetization in strongly nonequilibrium spin 
systems can be interpreted as magnon condensation, though the Holstein-Primakoff trans- 
formation for such nonequilibrium situations is not applicable. 

To exhibit condensation, unconserved quasiparticles require nonequilibrium conditions. 
Contrary to this, the condensate of conserved particles becomes depleted by nonequilibrium 
external perturbations. 
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Figure Caption 



Fig. 1. Qualitative phase diagram showing the sequence of states on the plane of the 
pumping field amplitude versus pumping time. 
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Figure 1: Qualitative phase diagram showing the sequence of states on the plane of the 
pumping field amplitude versus pumping time. 
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